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Abstract 

We present a unified picture of the interaction effects in dilute atomic quan- 
tum gases. We consider fermionic as well as bosonic gases and, in particular, 
discuss for both forms of statistics the fundamental differences between a gas 
with effectively repulsive and a gas with effectively attractive interatomic in- 
teractions, i.e. between a gas with either a positive or a negative scattering 
length. 
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I. INTRODUCTION 



One of the most important objectives of atomic physics in the last two decades, has been 
the achievement of Bose-Einstein condensation in a weakly interacting gas. It is, there- 
fore, not surprising that the observation of Bose-Einstein condensation in the magnetically 
trapped alkali vapors 87 Rb ffl], 7 Li M, and 23 Na ||, is considered to be one of the most 
exciting events of the last year. This is not in the least the case, because now that the 
successful road towards Bose-Einstein condensation has been found, one can safely assume 
that gases of atomic 85 Rb, 133 Cs, and 1 H will also be Bose condensed in the near future 
and that a large number of different experimental systems will be available to study various 
aspects of the condensation phenomenon in detail. 

Moreover, both lithium and hydrogen have stable fermionic isotopes which can be 
trapped and cooled in a similar manner as their bosonic counterparts. For atomic 6 Li 
this has already been achieved ||, but magnetically trapped deuterium ( 2 H) has not been 
observed yet because the loading of the trap cannot be accomplished in the same way as 
in the case of atomic hydrogen. This is presumably caused by the fact that i) deuterium 
binds more strongly to a superfluid helium film, ii) the surface recombination rate is much 
larger, and iii) the sample is contaminated with atomic hydrogen ||. Nevertheless, it is to 
be expected that these experimental difficulties can be overcome and that one will soon be 
able to study both these fermion gases in the degenerate regime. 

In view of the exciting developments mentioned above, it appears justified to present here 
an overview of the physical properties of all these weakly interacting quantum gases. We 
thereby intend to bring out most clearly the differences due to the Fermi or Bose statistics 
on the one hand, and due to the effective interatomic interaction being either repulsive or 
attractive on the other. This division in the focus of the paper is also reflected in its lay-out. 
In Sec. we present a general introduction into the physics of the atomic quantum gases 
of interest. In particular, we specify the approximations that are allowed for these quantum 
systems, and subsequently show how one can determine the effective interaction between 
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the atoms. Using the latter result, we consider in Sec. |(| the Fermi and in Sec. [IV] the 
Bose gases. Because we wish to bring out the different physical properties most clearly, we 
treat there only the homogeneous case. The experimentally more realistic case of a gas in a 
harmonic oscillator potential, however, does not lead to any qualitative differences, except 
for a Bose gas with an effectively attractive interaction. The latter possibility is therefore 
briefly discussed in Sec. |V|, where we also end the paper by pointing out some problems that 
remain for the future. 

II. INTERACTING QUANTUM GASES 

All atoms mentioned in the introduction are members of the first primary group of the 
Periodic Table of Elements and, therefore, have a magnetic moment fl that is equal to the 
sum of the magnetic moments of the electron spin s (with s = 1/2) and the nuclear spin 
i. Due to this magnetic moment, the atoms can easily be trapped by an inhomogeneous 
magnetic field and can be evaporatively cooled by a resonant microwave field that flips the 
magnetic moment. The successful experiments of last year have shown that this is a great 
advantage. However, the magnetic moment also gives rise to an important problem because 
two-body collisions in the gas can lead to a change in the component of the magnetic moment 
parallel to the magnetic field and thus to a decay of the atomic density in the magnetic trap. 

This decay mechanism seriously limits the lifetime of the gas, and can in principle be 
caused both directly by the magnetic dipole-dipole interaction V d and by the central (sin- 
glet/triplet) interaction V c in combination with the hyperfine interaction V h f = cihfS-i 
0. In general, the typical timescale for the former process is much larger than that for 
the latter, and the gas spontaneously spin-polarizes to a state in which the lifetime Ti ne \ is 
determined solely by the magnetic dipole-dipole interaction. There are roughly speaking two 
scenarios: At small magnetic fields B -C a/j///i e , where fi e denotes the electron magneton, 
the decay due to the central interaction can only be avoided by polarizing both the electron 
and the nuclear spin along the direction of the magnetic field, i.e. all the atoms are in the 
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spin state \m a = 1/2, mj = i). This corresponds to a so-called doubly spin-polarized atomic 
gas. However, for large magnetic fields B ^> ahf/fi e , the nuclear spin does not necessarily 
have to be polarized to obtain a relatively long lifetime and one can essentially trap any 
mixture of the states \m s = 1/2, mj) with an arbitrary projection of the nuclear spin. In 
this case we are dealing with an (electron) spin-polarized atomic gas. 

It is important to note that although the central interaction does not determine the 
lifetime of the (doubly) spin-polarized gas, it does determine the typical timescale r e i for 
elastic collisions. (See, however, Sec. |TT| for one exception.) Because r e i <C r ine i for not 
too low temperatures, we can consider the spatial degrees of freedom of the gas to be 
in equilibrium, whereas the spin degrees of freedom clearly are not. For timescales short 
compared to the lifetime, we are thus allowed to apply equilibrium statistical mechanics to 
the gas and our main problem is the accurate treatment of the elastic part of the central 
interaction, which for the (doubly) spin-polarized gases of interest is the triplet interaction 
V T . 

The determination of the effect of the interaction on the thermodynamic properties 
of a quantum gas is considerably simplified by the presence of two small parameters in 
the problem. In terms of density n of the gas, the range of the interaction ry, and the 
thermal de Broglie wavelength A = J^Kf^Jrnk^T of an atom with mass m, the two small 
parameters are the gas parameter nr v and the quantum parameter ry/A. Physically, the 
condition nr v <C 1 expresses the fact that we are dealing with a dilute system in which it 
is highly improbable that three particles are so close together that they can interact with 
each other simultaneously. As a result we can neglect three-body processes and only have to 
account for all two-body processes. Furthermore, the condition ry/A 1 shows that if two 
atoms interact, their relative angular momentum can only be equal to zero. In combination, 
we therefore conclude that for an accurate description of a quantum gas we only need to 
consider all possible two-body s-wave scattering processes. 
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A. Two-body transition matrix 



To explain how this may be achieved, we first of all introduce the appropriate 'grand- 
canonical' hamiltonian of the (doubly) spin-polarized atomic gas. Using the language of 
second quantization ]7|], it reads 

^ = E E (%a - Va)al a a^ a + ^EE V Q a l+?,a a l'^,a' a k',a' a k,a = H + H int , (1) 
a k aa ' kk'q 

where at Q creates and a^ a annihilates an atom in the (electron and nuclear) spin state |a) 
and the momentum eigenstate \k) of a cubic box with volume V and periodic boundary con- 
ditions. As always, these creation and annihilation operators obey the (anti) commutation 
relations [a^ a , at = 5 %z ,5 a>al and [a^ a> a Rl ^} ± = \a\^a\,J± = if the atoms are 
fermions or bosons, respectively. Moreover, V^= J dx VT{x)e~' bq x denotes the Fourier trans- 
form of the triplet interaction and = h 2 k 2 /2m + e a = + e a is the energy of the single 
particle states \k, a). Note that both the spin state |a) and its energy e a depend on the 
applied magnetic field B due to the Zeeman interaction V z = —pi - B. This dependence can 
for any particular atom easily be determined by a diagonalization of the spin part of the 
atomic hamiltonian, i.e of V h f + V z , but is not explicitly calculated here since it plays only a 
minor role in the following. Note also that in agreement with our previous remarks, we have 
only assumed that the spatial degrees of freedom are in equilibrium and have, therefore, 
introduced a different chemical potential fi a for each spin state \a). 

Next we consider the effect of the interaction H int on an 'initial' state \i) of two atoms. 
In particular, we take \i) to be the state in which the occupation numbers of all the single 
particle states are zero except for Ng, 2+ ^ a and Ng, 2 _£ a ,, which are equal to one. In lowest 
order in the interaction, the energy of this state is Ei = £g/2+ka + e K/2-ka' ^ we use ^hat 
the chemical potentials for this two-body problem are equal to zero. Applying standard 
perturbation theory, we find that this energy is shifted due to the interaction by an amount 

AE t = (i\H int \i) + J2 (i\Hint\f) B 1 B (f\H int \i) + . . . , (2) 
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where the 'final' state |/) is an eigenstate of H with energy Ef. Moreover, due to the 
possibility of scattering to other momentum states, the state \i) also acquires a finite lifetime 
Tj that up to second order perturbation theory is found from Fermi's Golden Rule 

1 2tt 

- = tE 6(Ek - E f )\(f\H M \i)\* . (3) 
Tj n j. 

Introducing the usual notation l/(x + iO) for the limiting procedure lim^o l/(x + ir]), both 
these results can be conveniently combined into 

AEl ~ Ir" = {ilHintli) + ^ ^ H ^ E ,-E f + tO {flHmtlt) + • • • • (4) 
The relevant matrixelements of the interaction Hi nt are easily evaluated and we obtain that 

ifi 1 | 1 1 

AB " " 2^ = V \ Vg + V f y "R 2 (fc2_fc« )/m + i0 V f'-« + ' ' ■ 

T + ^ E ^* a , ( g,_^ )/m + <0 ^ + ...)}, (5) 

which in the continuum limit V — > oo becomes 

- — = i J y(o) + / v(£ - £')— ^ — J w£' -£) + ... 

t w + / 7^ - feO r^s r L — - fc) + ■■■)} • (6) 

y J (27r) d h\k 2 - k' 2 )/m + lO )) 

The higher order contributions can be calculated in a similar way and we actually find 
that the complete perturbation series can be summed by introducing the so-called two-body 
T(ransition) matrix T 2B (k', k; E), which obeys the famous Lippman-Schwinger equation || 

T 2B (k>, k; E) = V(k' - k) + / ^ V(k' - k") -J- — T™{k\ k; E) (7) 

J E - h k" 2 /m + lO 

and is proportional to the amplitude for two atoms to collide with a kinetic energy E in 
the center-of-mass frame and change their relative momentum from fik to Hk'. Indeed, a 
comparison with Eq. (|6|) shows that the exact result is 

AEi - — = i {T 2B (k, k; h 2 P/m) T 8 aa ,T 2B (-k, k; h 2 k 2 /m)\ . (8) 



We therefore conclude that for a system of two atoms we can account for all two-body 
processes by first replacing in the hamiltonian of Eq. (|I|) the potential V$ by the transi- 
tion matrix T 2B ((k — k')/2 + q,(k — fc')/2; h 2 (k — k') 2 /4m) and then using only first-order 
perturbation theory to avoid a double counting of the effects of the interaction. 



B. Many-body transition matrix 

In a gas we are of course not dealing with only two atoms, but with a large number 
N = nV of atoms at a temperature T = l/ksft. This implies that we should repeat our 
previous calculation for a grand canonical ensemble of 'initial' states, in which the single 
particle states \k,a>) have the average occupation numbers N% a = (exp/3(eg a — /x Q ) ± 
and the chemical potentials determined from 

^ = ££%^£iv a . (9) 

The average zeroth-order energy in this ensemble is equal to 

^ = ££ (% a -^a)N la (10) 

01 k 

and its average energy shift and lifetime due to the interaction H int can now be calculated 
perturbatively by averaging Eq. (Q). Performing this thermal average is somewhat cumber- 
some, but ultimately we find 

ifi 1 

^■Ei — — = £ £ ^K/2+k,a^K/2-k,a' 

1 aa ' kk 

1 V P ^(k 2 - k' 2 )/m + *0 k '~ k 

neglecting between the curly brackets terms that are either of order 0(Vg (nry) 1 ^ 3 ) or of order 
0(V$ ry/A). (Diagrammatically speaking, we neglect the four 'bubble' diagrams compared 
to the two 'ladder' diagrams.) 



The appearance of the factor 1 =p Ng, 2+ % a =F a' * n ^his expression may appear 

surprising at first instance. It has, however, a clear physical meaning that becomes evident 
when we consider the contribution to the average lifetime from a scattering process with 
a 7^ a'. Such a contribution is given by 

2tt 1 „ Jh 2 k 2 _ h 2 k' 2 \ 2 / x 

^^/Z^ d l m m J \ V k-k'\ ^K/2+k,a^K/2-k,a' l, 1 T ^ K /2+k' ,a T ^ K /2-k> ,a> ) 

and consists of the sum of a loss rate 

E 5 - ^J-j l^-fe'| 2iV K/2+ifc ) a iV K/2-fe lQ ' (l T Nx /2+Pta ) (l =F Ng,^^) 

due to the collision of two particles, and a production rate 
2tt1 ^ J % 2 k 2 h 2 k' 2 \ |T/ 2 

^ y L " I m ml ' K/2+k,a^K/2-k,a'^K/2+k\a^K/2-k',a' 

due to the collision of two holes. For a scattering process with a = a' we can obtain a similar 
interpretation, although it is slightly more complicated in this case due to the interference 
effects. 

Following our discussion of the two-atom problem in Sec. [II A| , we see that we can 
now include all two-body processes by introducing a many-body T(ransition) matrix 
T^(P, k, K] E) that obeys 

T^(k',k,K ] E) = V(k'-k) 

J {^r E — h k" 2 /m + tO 

in the continuum limit. Using Eq. (0), this is equivalent to 

T^(k',k,K;E)=T™(k',h,E) 

T / * ^ ^/^^-^ Ig?^, t, *; £") , (13) 

■/ (27TJ £ - h z k" 2 /m + lO 

which is more convenient for our purposes since it offers the opportunity to explore the 
consequences of the small parameters nry and ry/A. Indeed, as a result of these small 
parameters the relevant momenta and energies are always small compared to K/ry and 
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h 2 /mry, respectively, and we can safely neglect the momentum and energy dependence of 
the two-body T matrix. Consequently, Eq. (|13|) is immediately solved by 

^*' giE)= i^S£(^ ' (14) 

where the quantity T 2B (0, 0; 0) = 4irah 2 /m is proportional to the (triplet) s-wave scattering 
length a and E aj0l i (K; E) is defined by 

Moreover, using this solution we find that in a good approximation 

A ^.-| = ^EE ^e/2 + M^?/ 2 - fe >(l T 6a*)T%?@, 0, i?; ft 2 ^/™) , (16) 

which can again be reproduced by first-order perturbation theory if we replace V$ in Eq. ([!]) 
by T a M ^(0, 0, k + fc'; fr 2 (£ - k') 2 /Am). 

With the last conclusion we have basically achieved our objective of summing all two- 
body s-wave scattering processes. In summary, we have shown that this requires the in- 
troduction of a many-body T matrix, which physically also incorporates the fact that the 
collisions of interest take place in a (gaseous) medium and are, therefore, influenced by the 
usual Fermi-blocking or Bose-enhancement factors [0. As a result, the transition matrix 
T^f l ^(k',k,K;E) acquires a dependence on the center-of-mass momentum TiK and the en- 
ergy E that is explicitly displayed in Eqs. (Ill) and (|T5|). This turns out to be important 
for the rest of the paper, where we use the many-body T matrix to arrive at an accurate 
(selfconsistent) theory of interacting quantum gases. 

III. FERMI GASES 

We start our discussion of the interacting quantum gases with the Fermi case. We will 
only consider spin-polarized atomic gases in which two hyperfine states are almost equally 
populated. For convenience, these two hyperfine states will be denoted by | f) an d | !), 
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respectively. The reason why we do not discuss doubly spin-polarized gases, is that in 
such gases the atoms can only scatter via p-waves due to the Pauli principle. As a result 
the (thermalizing) elastic collisions are dominated by the long-range magnetic dipole-dipole 
interaction V d and the spatial degrees of freedom only equilibrate on a timescale that is 
comparable to the lifetime of the gas, i.e. r e \ ~ r ine i. Furthermore, if one is able to cool the 
system by other means, the gas basically behaves as an ideal Fermi gas on timescales short 
compared to the lifetime and interaction effects are negligible (cf. Eq. (|16|)). 



A. The case a > 

We have argued that to account for all two-body s-wave scattering processes we must 
use the interaction 

Hint = i E Tl{ B (0, 0, k + k'; h\k - k'f/Am) a\ +fi / %l _^a p ^ . (17) 

kk'q 

However, this is not completely true because this interaction induces on average also a change 
in the hamiltonian Hq, due to the fact that the occupation numbers Ng a = (4 Q a fca) do 
not vanish. This so-called mean-field correction is clearly given by 

AH o = E E 4,a(t7 £ E Tu mk + k';h\k-kr/*rn)N p \ (18) 
a k { «¥« k' ) 

and, denoting the density of atoms in the spin state \a) by n a , is well approximated by 

AH = EE ( E T 2B (0,0;0K,} , (19) 

since the dominant contributions to the integration over the momentum hk' come from the 
region near the Fermi momentum Tikp^i = ^2m([i a t — t a i) (or from the thermal momentum 
h/A if the gas is not degenerate) where the many-body T matrix is almost equal to the 
two-body T matrix. The change in H can, therefore, be absorbed in a suitable redefinition 
of the chemical potentials. Introducing 

t*' a = t*«- E T 2B (0,0;0K, , (20) 
10 



the correct effective hamiltonian for the spin-polarized Fermi gases finally becomes 

H eff = EE (% a ~ »'aH a %* ~ VT 2B (0,6; 0)n T n 4 
a k 

+ ^ E 0, k + k'; h 2 (k - kf/im) a^a^ap^ , (21) 

kk'q 

where we have subtracted the constant VT 2B (0, 0; 0)n-|-n^ to avoid a double counting of the 
effects of the interaction that have already been included in the redefinition of the chemical 
potentials. Notice that we have implicitly assumed that the substitution /i a — > n' a is carried 
out also in the occupation numbers of Eq. (|T5|), since this makes the calculation of the 
many-body T matrix self consistent. 

From this effective hamiltonian we conclude that the spin-polarized gases are so-called 
Fermi liquids H]. In particular, the elementary excitations (or quasiparticles) of the gas 
have a dispersion relation given by eg — fi' a and scatter of each other with an amplitude 
given by the many-body T matrix. As a result the density of atoms in the two spin states 
is determined by 

1 



V 



E % 



a ' 



(22) 



with the occupation numbers N^ a equal to the Fermi distribution function (e@ x + 1) 1 



evaluated at e 



k,a 



k,a 

fi' a . Moreover, the pressure of the gas is simply 

-,2 



V 



^ EE ln ( 1 + e-^'^) + = EP« + Pint ■ 



V 



in 



(23) 



For a positive scattering length a, the interaction between the atoms is effectively re- 
pulsive and the second term in the right-hand side is positive. Interestingly enough, this 
does not imply that the gas is always mechanically stable. Indeed, using the hydrodynamic 
equations dn a /dt = —V • J a and J2 a dJ a /dt = —Vp/m, we find that that the two sound 
modes in the k clS £1X6 described by the coupled wave equations 

' <9p T /n T dp int /n i 
dpint/n^ dp l /n l 



d 2 




1 


dt 2 




m 







■V 2 




. n l. 



(24) 
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We thus conclude that the velocities of the sound modes are equal to the square root of the 
eigenvalues of the matrix in the right- hand- side of Eq. fl2~4]), and that both these eigenvalues 
need to be positive for the sound velocities to be real and the gas to be mechanically stable. 

In the case of a nondegenerate Fermi gas we can apply Boltzmann statistics to obtain 
Pa — n a kBT. The condition of mechanical stability then reduces to n^n^a 6 < (a/A) 4 /4 
and requires rather low densities. In the more interesting case of a degenerate Fermi gas, 
however, we have 

1 / fW \ 2/3 
Pa = n a k B T- (3J-n a A 3 ) , (25) 



5 V V 2 

neglecting terms of order 0(i/(fcp jQ A) 4 ). The condition on the densities then becomes 
r^n^a 6 < 7r 2 /2304, which is much less restrictive. Note that in the mechanically unstable 
region, the gas is in general unstable against density fluctuations and will phase separate 
into a dense (liquid or solid) and a dilute (gas) phase. This process of phase separation is 
known as a spinodal decomposition and the density and temperature conditions at which it 
first occurs, lie in our case on a so-called spinodal surface. 

Apart from the first-order phase transition implied by the above instability, we do not 
expect any other phase transition to occur in the gas on the basis of our effective hamiltonian. 
In particular, we do not expect a phase transition due to quantum degeneracy, because for 
that we need an attractive interaction between the quasiparticles as will be explained shortly. 
In principle, however, it is possible that a positive s-wave scattering amplitude induces a 
negative p-wave scattering amplitude due to the Kohn-Luttinger effect ||10|| . In this manner 
a phase transition to a superfluid state might occur after all. Although it is interesting to 
analyze this possibility in more detail, we will not do so here since it is beyond the scope of 
the present paper. 



B. The case a < 



The prospects for achieving a phase transition in the degenerate regime are more promis- 
ing for spin-polarized Fermi gases with a negative scattering length, which actually applies 
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to both atomic deuterium 11] and atomic 6 Li |T2[ . Nevertheless, for this phase transition 



to be experimentally observable, it must take place in the (meta)stable region of the phase 
diagram where the gas is mechanically stable. Since the stability analysis of the previous 
section is insensitive to the sign of scattering length, this implies in the degenerate regime 
that the densities of the two spin states have to fulfill the requirement n^n^a 6 < 7r 2 /2304. 
In contrast to the case of a Fermi gas with a positive scattering length, the latter condition 
now does not rule out the possibility of the formation of a superffuid state at sufficiently 
low temperatures. 

The physical reason for this is that with effectively attractive interatomic interactions the 
gas can form a Bose condensate of Cooper pairs, in the same way as an electron gas can form 
a condensate of Cooper pairs in the BCS theory of superconductivity |J. Mathematically, 
the instability of the gas towards the creation of Cooper pairs is signalled by the fact that the 
scattering amplitude diverges for two quasiparticles at the Fermi energy = h 2 k 2 F /2m = 
(h 2 kp^/2m + h 2 kp J2m) /2 with opposite momenta and 'spin' [Tj|. Using Eqs. ([TJP and 
(|T5|) and neglecting the small imaginary part of 2(0; 2e F ), we thus find that the critical 
temperature of the BCS transition obeys 

1 m m r dk Nt(k) + N\(k) n 



T# B (0,0,0;2e F ) Airah 2 J (2vr) 3 2e F -h 2 k 2 /m 

where V denotes the Cauchy principal- value part of the integral. 

The properties of this so-called linearized BCS gap equation are well-known from the 
theory of superconductivity. First, this equation only has a solution if a < 0. Second, the 
highest critical temperature for a given total density n is obtained if nj = nj_ or equivalently 
if 5e F = \h 2 k 2 F ^/2m - % 2 k 2 F j2m\ = 0. This can be understood most easily by noting that 
for a nonzero value of 5€f there are less particles at the Fermi surface that can be paired, 
which suppresses the critical temperature. Third, if the densities in both spin states are 
equal, the gap equation can be solved analytically |IJ| and leads to 

^IH-s^- 1 } • (27) 
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For unequal densities an analytical treatment is not feasible, but a numerical solution of 
Eq. fl26|) shows that nonzero critical temperatures are only possible if the 'spin polarization' 



n\ —ni\/(ni+ni) is less than 3kBT c /2eF fl^l - Qualitatively, this is a result of the fact that 



if SeF is of the order of fceT c , it is no longer energetically favorable to form Cooper pairs 
and the gas is in a normal state even at zero temperature. 

From an experimental point of view, we thus conclude that the most favorable conditions 
for the observation of the BCS transition in a spin-polarized Fermi gas are obtained if n-f = 
ni = n/2. Under these conditions the mechanical stability of the gas at the temperatures of 
interest requires that n\a\ 3 < 7r/24 or equivalently that kp\a\ < 7r/2. For atomic deuterium 
the triplet scattering length is equal to a ~ —6.8 a (a is the Bohr radius). At realistic 
densities of say 10 12 cm -3 we thus find that kp\a\ ~ 1.1 x 10~ 3 <C tt/2 and the stability 
condition is easily fulfilled. However, from the exponential dependence in Eq. (^7|) we 
immediately see that the critical temperature is completely out of reach with present-day 
cooling techniques. Fortunately, this conclusion does not hold for atomic 6 Li because the 
triplet scattering length now has the anomalously large value of a ~ —4.6 x 10 3 a^. For the 
same density of atoms we thus find that kp\a\ ~ 0.74 < tt/2 and that the critical temperature 
is only 29 nK, which is rather close to the temperatures that recently have been obtained 
with the bosonic isotope 7 Li H . Therefore, it appears that 6 Li is a very promising candidate 
for the observation of the BCS transition in a weakly interacting gas. 

Because of this possible application, we also want to mention briefly how we can arrive 
at an accurate description of the gas below the critical temperature. Although it is not 
difficult to treat the general situation with ^ we concentrate here on the most 
important case of an equal population of the two hyperfine states. Moreover, we also only 
discuss the derivation of the elementary excitations of the superfluid state (which is sufficient 
for the most interesting equilibrium properties of the gas) and do not consider the effective 
interaction between these quasiparticles. Introducing the convenient notation // = fi' a — e a 
for the renormalized chemical potential (or Fermi energy e^) of both spin states, our task is 
therefore reduced to including into the free hamiltonian of the quasiparticles, i.e. into 
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H e ff -EE - A*04, a a^ - VT2B & o; °)t ' (28) 



a fe 

the effect of the Cooper pairs. 

To achieve this, we must realize that BCS theory is roughly speaking the theory of a 
Bose-Einstein condensation of Cooper pairs. Therefore, the order parameter of the phase 
transition is the expectation value of the annihilation operator for a pair of atoms. (See 
Sec. [IV] for a treatment of Bose-Einstein condensation in an atomic Bose gas.) As a result, we 
do not only have nonvanishing average occupation numbers (o^ a d^ a ) but also nonvanishing 
values for (o-^ a a _j: a i)- Note that in our case we are dealing with a pairing between atoms 
with different spin states and the latter average is only nonzero for a ^ a'. To see what 
the influence of the anomalous averages (a^^a_^,) is, we briefly return to the hamiltonian 
of Eq. ([!]) . Taking for the sake of clarity = Vg, the lowest order average correction on the 
hamiltonian Ho due to the interaction now becomes (cf. Eq. (|I9"D) 

(29) 



AH = EE a l a %a E V n «' \ + E { A o4t a -fc| + A S a -fc,| a fc,t} 

a k Uv<* J k 



if we define the BCS gap parameter A by 

A o-^E V o( a -k,i%J ■ ( 3 °) 

k 

From our previous discussion of the normal state of the gas we know that we can to a good 
approximation account for all two-body processes by replacing in Eq. ( U% the potential 
Vq by T 2B (0, 0;0). It is important, however, that we do not perform this substitution in 
Eq. ( |30|) because the BCS theory is already going to sum all the other relevant two-body 



processes for us [15 



We thus arrive at the conclusion that the free hamiltonian of the quasiparticles, including 
the effects of a condensate of Cooper pairs, is given by 

H o" = E E ( e A? - ^>{ a a k,a + E { A o4, T a -fc,l + A *o a -k,l a kA} 

a k k 

-V l -^-VT™(0A0)^, (31) 
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where we have again subtracted a constant term to avoid a double counting of the interaction 
effects. Being quadratic in the creation and annihilation operators, this hamiltonian can be 
diagonalized by means of a Bogoliubov transformation and we find that the Bogoliubov 



quasiparticles have an energy dispersion hu% = y(e^ — /x') 2 + |A | 2 which for nj = is 
independent of their 'spin'. Moreover, after this diagonalisation the average density n and 
order parameter A are easily calculated and result in the equation of state 

and the BCS gap equation 

1 1 _ 1 - N % t - N % , 

v. + v^ ik - = °> < 33 > 

respectively. Here, the average occupation numbers of the quasiparticles are again denoted 
by N% , but they are now equal to the Fermi distribution function evaluated at hu^ and in 
fact independent of a. 

The equation of state and the BCS gap equation are, at a given density and temperature, 
two equations for the two unknown quantities fj,' and |A |. However, before we can actually 
solve these equations we first have to resolve the ultraviolet divergence in the gap equation. 
This divergence is a consequence of the fact that we have neglected the momentum depen- 
dence of the potential and used V$ = Vg. Making use of the same approximation in the 
Lippman-Schwinger equation Eq. (^) and of the fact that the Fermi energy \il = Ti 2 k 2 F /2m 
is always much smaller than h 2 /mry, we find that the two-body T matrix obeys 

1 1 11^1 

— — ~ — = 1 V (34) 

T 2B (0,0;0) T 2S (0,0;2/i') V 8 2{e % - //) 

in this case. We thus observe that the divergence in the BCS gap equation can be cancelled 
by a renormalization of 1/Vg to 1/T 2B (0, 0; 0) = m/Anah 2 . In this manner we then obtain 
the following gap equation 

Anah 2 + V\ { 2%u % ' 2^-/)/"°' (35) 
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which is free of divergences. Furthermore, taking the limit Ao — > to obtain an equation 
for the critical temperature T c , we exactly reproduce (after taking also the thermodynamic 
limit and substituting jj = ep) the result of Eq. (|26| ) which shows the consistency of our 
approach and, in particular, also that the BCS theory in essence indeed sums all two-body 
scattering processes. 



IV. BOSE GASES 

We now turn our attention to the atomic Bose gases, which have caused so much excite- 
ment in the atomic physics community after the experimental observation of Bose-Einstein 
condensation in atomic 87 Rb, 7 Li and 23 Na vapours. As in these experiments, we consider 
here only the doubly spin-polarized case. This implies that there is only one spin state in 
the problem and we will, therefore, in this section completely suppress the spin indices in 
our notation. 



A. The case a > 



Again, we first consider a gas with an effectively repulsive triplet interaction, which 
applies for example to atomic hydrogen |§, atomic 87 Rb [|l6j and atomic 23 Na |L7]. The 



appropriate effective hamiltonian for the doubly spin-polarized atomic Bose gases is, in 



analogy with the discussion in Sec. Ill A, given by 



H eff = E (<* - ^>h - VT 2B & 0; 0)n 2 

k 

E tMB & ^+ - £') 2 /4m) ' ( 36 ) 

kk'q 

with a renormalized chemical potential \J equal to 

y! = ^-2nT 2B (0,0;0) . (37) 
Looking at this hamiltonian we immediately see that the density of atoms is determined by 
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" = 75>s (38) 

k 

where the average occupation numbers N? are now equal to the Bose distribution function 
(e^ x — evaluated at eg — //. In addition, the pressure of the gas is 



V 



^Vlnfl- c -«^>) + , (39) 

k 



which for a gas with a positive scattering length does not show a mechanical instability, 
i.e. dp/dn is always greater than zero for densities such that nry <C 1 and our theory is 
applicable. 

On the basis of Eq. fl38|) we thus expect the gas to Bose condense at the same critical 
temperature 

^ [cm) (40) 

as the ideal Bose gas. For an interacting Bose gas with positive scattering length the order 
parameter corresponding to this phase transition is the expectation value of the annihi- 



lation operator for an atom QTSj| . If we want to describe the gas also below the critical 
temperature, we therefore need to consider the influence of a nonvanishing average (aj) on 
the effective hamiltonian. Focussing again on the derivation of the elementary excitations 
in the superfluid state, this is most easily achieved by first carrying out the substitution 
a k ~ > ^k o( a o) + a k m t ne hamiltonian H e ^ and then keeping only all terms that are of first 
or second order in the creation and annihilation operators. In this manner we find in a good 
approximation that 

# e// = -V (V 2B (0, 0; 0)(n') 2 + T MB (0, 0, 0; 0)n'n o + ^T MB (0, 0, 0; 0)n 2 



+ vVn (-// + n T MB (0, 0, 0; 0)) (a 5 + at) + £ (eg - //)ata g 

k 

+ £ { (2n T MB (0, 0, 0; 0)) a\ H + ^n T MB (6, 0, 0; 0) (aja^ + a_ m ) } , (41) 

k 

if we choose (a^) to be real and introduce the condensate density n = (a^) 2 /V and the 
density of noncondensed atoms n' = n — hq. The most difficult part of the calculation is 



again to avoid double countings of the interaction effects. Besides the subtractions in the 
first term in the right-hand-side, this also requires that the relation between \i and /jf is 
changed into fi = 2n'T 2B (0, 0; 0) + //, which reduces to Eq. ( p7|) if the condensate density 
vanishes. 

To proceed, we first note that if we want (aj) to be the total expectation value of the 
original operator ag, we must require that the expectation value of the new operator 
(i.e. after the substitution) disappears. This is achieved by putting fjf = n T MB (0,0,0; 0), 
because then the terms in the hamiltonian linear in the creation and annihilation operators 
vanish. As a result, the chemical potential is determined by 

fi = 2(n - n )T 2B (0, 0; 0) + n T MB (6, 0, 0; 0) , (42) 

which for the approximation that we are using is just the famous Hugenholtz-Pines relation 
fT9|| . Because our theory obeys this relation, we expect that the dispersion relation hu^ of 
the quasiparticles is linear at long wavelengths. This indeed turns out to be the case, since 
a diagonalization of the above hamiltonian by means of a Bogoliubov transformation shows 
that hu % = y / e| + 2n T MS (0,0, 0; 0)e^ and that 

Ho" = E 

k 

-V (t 2B (6, 0; 0)(n') 2 + T MB (6, 0, 0; 0)n'n + ^T MB (6, 0, 0; 0)ng) , (43) 

where fot and denote the creation and annihilation operators for the Bogoliubov quasi- 
particles and we have neglecting between the brackets additional terms that are a factor of 
order 0(V% a3 ) smaller ||18|| . Moreover, using this Bogoliubov transformation we can now 
obtain also for the superfluid state the pressure 

P = -^£ ln(l-e-^) 

k 

+T 2B (6, 0; 0)(n') 2 + T MB (6, 0, 0; 0)n'n + ^T MB (6, 0, 0; 0)n 2 (44) 
and the equation of state 
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„ „ y [ eg + «oT^(0,0,0;0) eg + n T^(0, 0, 0; 0) - hu £ \ 
° + 1 I K 2 ^ J ' 

which determines the condensate density at a fixed density and temperature. Note that in 
the equation of state the average occupation numbers of the quasiparticles N% are equal to 
the Bose distribution function evaluated at Twor- 

This almost completes our discussion of the interacting Bose gas with a positive scattering 
length. There, however, remains one point that we need to address, namely our theory is 
not selfconsistent yet. This is because we have not yet mentioned how the many-body T 
matrix is determined below the critical temperature T c . In general, this is a complicated 



problem due to the presence of infrared divergences in the theory of a dilute Bose gas |2T[ . At 
very low temperatures such that naA 2 3> 1 the treatment of these divergences requires more 
advanced methods then the ones we are using here. Therefore, we consider only the opposite 
limit naA 2 <C 1, which is in fact the most relevant one for experiments. In this regime the 
average kinetic energy ksT of the atoms is much larger than the average interaction energy 
Ananh 2 /m and the Bogoliubov dispersion is for all practical purposes well approximated by 
huj: ~ + n T MB (0, 0, 0; 0) = + In Eq. this means that we can neglect the 
anomalous terms proportional to cb\a}_^ or a_^a^.. Hence, we now basically deal with the 
hamiltonian Yl% ( e £ + ^') a \ a k anc ^ the many-body T matrix is again just given by Eqs. (|HD 
and flTSp. Finally, we have to deal with one last subtlety. In the way we have defined the 
many-body T matrix, it describes the scattering of quasiparticles. However, we are dealing 
with a Bose condensate of particles and we actually need their scattering amplitude. In 
our case this, fortunately, only leads to a slight shift E — > E — 2/x' in the energy of the 
many-body T matrix. We thus arrive at a fully selfconsistent theory of the superfluid phase 
in the regime naA 2 -C 1, if in all formulas below Eq. (flUD, T MB (0, 0, 0; 0) is replaced by the 
many-body T matrix T MB (0, 0, 0; -2fi'). 

The transition matrix T MB (0, 0,0; — 2u') is almost equal to T 2B (0, 0; 0) = 4irah 2 /m for 
most temperatures below T c , but decreases to zero for temperatures near the critical tem- 
perature due to the fact that 2(0; —2//) diverges in the limit u' — > 0. Therefore, our self- 
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consistent approach reduces to the conventional Bogoliubov (or Popov) theory of a weakly 
interacting Bose gas for most temperatures. However, near the critical temperature it re- 
solves a long-standing problem of the Bogoliubov theory, which actually predicts a first-order 
phase transition |22| instead of the second-order phase transition expected from the theory 



of critical phenomena. Of course, using the many-body T matrix instead of the two-body 
T matrix does not include all the critical fluctuations that are of importance near the crit- 
ical temperature. Nevertheless, using the many-body T matrix we incorporate, at least 
qualitatively correctly, the fact (known from renormalization group theory ||23|| ) that the 
interaction between the condensate particles should vanish exactly at the critical tempera- 
ture. It is rather interesting that this can already be achieved by including only all two-body 
scattering processes. 

B. The case a < 

We now turn to the physics of a Bose gas with effectively attractive interactions, which 
is realized for instance in atomic 7 Li |24| and atomic 133 Cs ||25|| . This is an interesting topic 
because a thoughtless application of the theory of Sec. |1V A| would result in a dispersion of the 



Bogoliubov quasiparticles that is equal to Tiuj^ = y e| — (87r\a\n h /m)e^ and therefore purely 
imaginary at long wavelengths. From this simple argument we can thus already conclude 
that Bose-Einstein condensation in the canonical sence, i.e. with the order parameter (a^), 



will not occur if the scattering length is negative [26|. However, from the experience with 



fermionic gases gained in Sec. HI Bj , we might expect a BCS-like phase transition to take 



place instead. This would be signalled in Eq. (|3q) by a divergence of the scattering amplitude 
of two quasiparticles at energy // and with opposite momenta, or by (cf. Eq. (|26|)) 

J„ /Jg- *® =0 (46) 

T A/B (0,0,0;2/i') T 2B (0,0;0) J (2vr) 3 ^-h 2 k 2 /2m 

with of course N(k) the Bose distribution function evaluated at h 2 k 2 /2m — jj! . Analyzing 
this condition, we find that the phase transition indeed occurs at a temperature Tbcs = 
T (l +0(\a\/ A )) which is only slightly above the critical temperature of the ideal Bose gas. 
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Following the same procedure as in Sec. |I 1B|, we can easily show that below this critical 
temperature the equation of state becomes 

V 2 ? \ hut k+ 2hu- f 1 j 



k 



k k 



and that the BCS gap equation is well approximated by 



TTl 1 N? 

+ 77 £ d: = o. (48) 



Aiiah V ^ huj 

k 

Here the energy of the Bogoliubov quasiparticles is hu^ = <J (eg — fi') 2 — |Ao| 2 and their 
average occupation numbers iVg are equal to the Bose distribution function evaluated at this 
energy. These equations again determine // and |A | at a fixed density and temperature. 

Notice that in contrast with the fermion case, the dispersion huj? has a minus sign in 
front of |A | 2 . This has important consequences because it implies that if we lower the 
temperature at a fixed density, both |A | 2 and fi' increase in such a manner that the gap 
in the Bogoliubov dispersion decreases. Hence, at a second temperature T BEC < T B cs 
the gap closes and the number of particles in the zero-momentum state diverges, which 
signals a Bose-Einstein condensation. Below that second critical temperature we always 



have |A | = — y! and the dispersion hu^ becomes equal to «/e| — 2/i'eg which, interestingly 
enough, is exactly the same Bogoliubov dispersion as for a Bose-condensed gas with positive 
scattering length. Moreover, Eqs. ([|7]) and ( f48D turn into 

and 

m rin 1 ^ Nr 

+ = < 50 > 

determining now // and the condensate density n . It is important to mention here that the 
possibility of a Bose-Einstein condensation in a gas with effectively attractive interactions, 
in the way that we have just seen, does not contradict the argument against Bose-Einstein 
condensation presented in the beginning of this section. The reason is that in that argument 
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we assumed that a Bose-Einstein condensation would be associated with the order parameter 
(ag) and did not consider the possibility that it would be associated with the condensate 
density no, as in the ideal Bose gas pf] . 

It thus appears that a Bose gas with a negative scattering length offers the exciting 
possibility of observing both a BCS transition and a Bose-Einstein condensation. However, 
we have up to now not considered the mechanical stability of the gas in these two superfluid 
phases. This is clearly of the utmost importance, because Eq. (|39D shows that for a negative 
scattering length dp/dn does not always have to be positive. Indeed, calculating the spinodal 
line from dp/dn = 0, we actually find the unfortunate result that the gas is always unstable 
at T BC s- We are thus forced to conclude that the two second-order phase transitions that 
we have found above are always preempted by a first-order gas-liquid or gas-solid transition 
and are, therefore, in a homogeneous system unobservable. 

V. DISCUSSION AND OUTLOOK 

In this paper we have presented a unified picture of weakly interacting atomic gases, 
focussing on the differences that arise due to the statistics of the atoms and due to the 
interatomic interaction being either effectively repulsive or attractive. However, we have 
restricted our discussion to the homogeneous case, whereas the experiments are up till now 
always performed in an inhomogeneous magnetic trap. Fortunately, the experiments are also 
always performed at temperatures such that ksT ^> hu, where hu is the energy splitting 
of the (harmonic oscillator) trapping potential. Under these conditions we are generally 
justified in using a local-density approximation, which means that we consider the gas as 
homogeneous in each point in space. In this sense the theory presented above can also 
be applied to an inhomogeneous situation and we, in first instance, do not anticipate any 
qualitative differences. Nevertheless, for quantitative predictions of for example the critical 
temperature and the density profile in the trap, it is of course essential to include the effect 
of the trapping potential. 



23 



Moreover, for a Bose gas with negative scattering length there is also a qualitative dif- 
ference. In Sec. [IV B| we argued that for a homogeneous gas a Bose condensate is always 



unstable in this case. However, first Hulet p7| and subsequently Ruprecht et al. [BHJ have 



argued that in an inhomogeneous situation one can have a metastable condensate if the 
number of particles is sufficiently small. This was recently put on a more firm theoretical 
basis by Kagan et al. |29| and by one of us [ 30| . In particular, the last author has shown 



by a calculation of the decay rate of the condensate due to thermal fluctuations, that a 
long-lived metastable condensate can exist for temperatures sufficiently close to the critical 
temperature. Although this result appears to explain the recent observation of Bose-Einstein 
condensation in atomic 7 Li, our theoretical understanding of these experiments is still in- 
complete. The most important problem that remains to be solved is the stability of the gas 
just above the critical temperature. Clearly, on the basis of a local-density approximation 
we would conclude, as in Sec. [IV B| , that the gas is unstable and that also in an inhomo- 
geneous situation Bose-Einstein condensation is unobservable. However, the local-density 
approximation is invalid near the critical temperature due to the large value of the (homo- 
geneous) correlation length. A thorough discussion of this point, therefore, requires a theory 
that goes beyond this approximation and is evidently of great interest in view of the Rice 
experiments with 7 Li ||. 

Also in the case of a Bose gas with positive scattering length, interesting theoretical 
problems remain. One of them is the correct treatment of the infrared divergences in the 
regime naA 2 3> 1, which, due to the effectiveness of evaporative cooling techniques, is also 
anticipated to be experimentally accessible. As already mentioned, we presumably need more 
advanced methods than discussed here to accurately perform a resummation of the divergent 
contributions to the thermodynamic quantities of interest. Indeed, we are presently carrying 
out a renormalization group study of the interacting Bose gas in this regime, because by 
this method infrared divergences are automatically avoided. Moreover, a resolution of these 
divergences is also of importance for the theory of two-dimensional atomic Bose gases, such 
as for instance doubly spin-polarized atomic hydrogen adsorbed on superfluid 4 He films. The 
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reason for this is that in a two-dimensional geometry, the divergences are so important that 
they actually prevent the possibility of a Bose condensate (in the thermodynamic limit). 
Nevertheless, the gas is expected to undergo a phase transition to a superffuid state at 
sufficiently low temperatures. In our opinion, formulating an accurate microscopic theory 
ST| for this so-called Kosterlitz-Thouless transition is an important challenge for the future. 

Finally, we would also like to mention the possibility of studying mixtures of fermionic 
and bosonic atomic gases. One exciting possibility is a mixture of doubly spin-polarized 
atomic 6 Li and 7 Li. In this manner one would create the dilute analogue of 3 He and 4 He 
mixtures. Such a mixture of doubly spin-polarized atomic 6 Li and 7 Li is rather interesting 
because, as we have seen in Sec. |HT[ the 6 Li atoms essentially do not interact with each 
other but only with the 7 Li atoms. However, due to the latter interaction they can now 
also be evaporatively cooled on a timescale which is much smaller than their lifetime. It 
therefore seems worthwhile to extend the methods discussed in this paper to this system 
and in particular study Bose-Einstein condensation in the presence of an essentially ideal 
Fermi gas. 

It is clear from these examples, and from many others that we unfortunately are not able 
to discuss here, that we are heading towards an exciting future, in which no doubt many 
surprises will await us. We hope that this paper will be of some use to people that intend 
to contribute to this rapidly developing area of atomic physics. 
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